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I. INTRODUCTION 

The question of the screening of nuclear pairing is a 
long standing problem. In the past it has mostly been 
considered in neutron matter (see Ref. P| for a list of 
references). However, quite recently also conjectures 
about the possibility of anti-screening, i.e., enhancement 
of pairing in finite nuclei due to surface vibrations have 
been put forward In Ref. 0 we also contributed to 
the question of the screening of the bare nucleon-nucleon 
pairing force in pure neutron matter and found consis¬ 
tently with the results of other authors that this polar¬ 
ization quenches quite strongly the neutron pairing in the 
S' = 0, T = I channel. However, as far as we know, simi¬ 
lar investigations have never been carried out in nuclear 
matter. Being interesting in its own right, our interest 
here is to try to make a qualitative link to hnite nuclei. 

Despite the fact that the validity of the local den¬ 
sity approximation (LDA) may be questioned due to 
the large coherent length of the Cooper pairs, it still 
yields, for example, for the pair correlation energy, semi- 
quantitative agreement with microscopic Hartree-Fock- 
Bogolyubov (HFB) calculations in hnite nuclei Q. Let 
us therefore resume our present knowledge about nuclear 
pairing in different channels: as already mentioned, prac¬ 
tically all calculations in the spin singlet channel of neu¬ 
tron matter indicate that polarization quenches pairing. 
What about ^ = 0, T = 1 pairing in hnite nuclei? In 
an important very recent contribution to this subject it 
has been shown that the bare nucleon-nucleon (NN) 
force yields already ^ 50% of the measured gap in the 
tin isotopes. It must be mentioned that in that calcu¬ 
lation the usual effective density dependent A:-mass with 
m*/m ~ 0.7 at saturation has been employed. This is in 
line with the common lowest order approach to pairing 
in nuclear physics where for example with the successful 
Gogny DIS force |2|in HFB calculations implicitly also 
the fc-mass is used Q- It further should be mentioned 
that within the same approach and specihcally with the 
Gogny DIS force the S = Q, T = 1 gap in symmetric 
nuclear matter is very close to the one obtained with 
the bare force for densities p < po/5 whereas it drops 
off quite a bit slower for densities p > po/S- At satu¬ 


ration (p = po), the bare force yields almost negligible 
gap whereas DIS still yields A ~ 0.5 MeV (see Fig. 2 in 
Ref. 0 )- In a LDA picture it is therefore not unreason¬ 
able that with the bare force the gap is about a factor 2 
smaller than the one from this DIS-force |^. The DIS 
force has therefore pairing properties which are not very 
far from the bare one. Remember that the weak coupling 
formula predicts an exponential dependence of the gap on 
the force (which however may be somewhat questioned 
in finite nuclei due to the discreteness of the spectrum). 
Indeed it has been shown recently that a renormalized 
bare force (riiow-k) which is phase shift equivalent for low 
energies (< 300 MeV) has shape and magnitude very 
similar to the Gogny DIS force Q- The fact that DIS 
is close to the bare force in the S = 0, T = 1 pairing 
channel has also been noticed by Bertsch et al.|i| in their 
investigation of the neutron halo in ^^Li. 

Concluding this discussion, we can expect that medium 
renormalization of the bare S' = 0, T = 1 pairing inter¬ 
action should yield some additional attraction in nuclear 
matter. This is indeed what is claimed to be the case 
in Ref. M and it is also what we will hnd in present 
study of infinite nuclear matter, employing exactly the 
same approach as the one used before in neutron matter 
0. However, as we will see later, the effect is dramati¬ 
cally strong. In finite nuclei, for instance in the presently 
very actively studied N ~ Z nuclei, we have the fur¬ 
ther very interesting neutron-proton [np) S = 1, T = 0 
pairing channel. It is generally believed that the pairing 
force in this channel should be of similar strength, may 
be a little stronger, than the one in the S = 0, T = 1 
channel M- In Refs. 0111 in the gap in the S' = I, 
T = 0 channel, calculated with the bare force (and k- 
mass) in inhnite symmetric nuclear matter, is given as a 
function of density, see e.g.. Fig. 1 of Bef. and Fig. 1 

of Ref. 0 ■ Because of the stronger attraction of the bare 
force in the deuteron channel the np gap turns out to be 
much larger than the neutron-neutron (nn) or proton- 
proton (pp) one. For example at maximum the np gap is 
about 8 MeV (!) whereas it is 2.5 MeV for the nn gap. 
Even at saturation still is of the order of 2 MeV 
(!). Glearly, the np gap in finite nuclei would turn out 
much too large, if the bare force was employed. Screening 
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therefore should give additional repulsion in the S' = 1, 
T = 0. We therefore have from the experimental side, 
and keeping in mind that LDA should at least give the 
right trend when going from the infinite homogeneous 
case to finite nuclei, definite predictions what the inclu¬ 
sion of polarization in a nuclear matter calculation should 
give, as a trend, in the S = 0, T = 1 and S = 1,T = 0 
channels: additional attraction (over the bare interac¬ 
tion) in the former and repulsion in the latter. We will 
see in how far these expectations are fulfilled by the cal¬ 
culations. 


II. BUBBLE SCREENING 



FIG. 1: Pairing interaction with screening in the RPA approx¬ 
imation. The dashed lines represent the Gogny interaction, 
the wiggly line the p-h residual interaction resummed to all 
orders. All vertices are to be understood as anti-symmetrised 
(not shown explicitly). 


The superfluid phase of a homogenous system of 
fermions is characterized by the pairing field Ai. (u;), 
which is the solution of the generalized gap eauation[l4L 

Hi in 


Ak{uj) = 'Y2 j 


( 1 ) 


where V is the sum of all irreducible NN interaction terms 
and Fk{uj) is the anomalous propagator. Dealing with a 
strongly correlated Fermi system one expects the medium 
corrections to play a crucial role. The gap equation by 
itself embodies the full class of particle-particle (p-p) lad¬ 
der diagrams just taking the bare interaction for V 0 
The same kind of correlations are incorporated in the 
propagators if the self-energy is approximated by the 
Brueckner-Hartree-Fock (BHF) mean field 0 From 
this side additional p-p corrections are unlikely to im¬ 
prove the predictions since the pairing rapidly vanishes 
at high density. On the contrary, particle-hole (p-h) cor¬ 
relations should play an important role at low density 
where the pairing gap is expected to take the largest 
value. But again p-h contributions have to be treated on 
equal footing, in the vertex corrections as well as in the 
self-energy in view of possible strong cancellations. This 
was the main concern of the study presented in Ref.Q 
for ^S'o nn pairing in neutron matter, further extended 
in a preliminary work to nuclear matter in 0 . The ex¬ 
pansion of the interaction block V and the self-energy 
S were both truncated to the second order in the inter¬ 
action. This approximation turns out to be reasonable 
for the self-energy, whereas it is only indicative for the 
trend of the screening effects due to the interaction. In 
fact it was shown that the screening on the nn pair¬ 
ing is quite different according to whether the medium is 
made out of neutrons or neutron plus protons. Quantita¬ 
tive predictions require that the full RPA bubble series, 
at least, should be summed up in order to calculate the 
screening interaction, being aware that even this approx¬ 
imation could be not enough to remove the singularity 
associated to the low-density instability of nuclear mat¬ 
ter. 

According to the preceding discussion, we extend the 
study of Ref. 13 including in V the full bubble series. It 


is convenient to develop the latter as shown in Fig. ^ 
where the block (c) is the correction to the results of 
Ref.Q- The splitting of the bubble series into the dia¬ 
grams (b) and (c) enable us to disjoin the consideration of 
the two external vertices connecting one particle line with 
one hole line from the full p-h vertex (wiggly line). So 
doing, the former two vertices can be calculated taking 
exactly into account their full momentum dependence, 
whereas the latter can be restricted to transitions around 
the Fermi momentum. In fact it is well known that the 
magnitude of the gap is quite sensitive to high momenta 
transitions, i.e. the short range part of the nuclear force 
0 , whereas p-h excitations are only important near the 
Fermi surface. 

The p-h bubble series is summed up by the Bethe- 
Salpeter equation 0 


^J^T^(9i>92,g) 


/ ^4 

xVs,Tiqi,q3,q)V^^^{q3,q2,q). ( 2 ) 


Here q = (w, q) and S, T denote the total spin and isospin 
in the p-h channel. q) = —iG{q 3 )G{q + ga) is the 

polarization insertion|l8l|. Since, as we said, only small 
energy-momentum transfers significantly contribute to 
the pairing interaction around the Fermi energy, it is 
likely a valid approximation to replace the p-h residual 
interaction with the effective interaction V expressed in 
terms of Landau parameters: 


V = f + 9 {o- ■ o-') + f {t ■ t') + g' (a- ■ cr') (r • r'). (3) 


The remarkable advantage of this approximation is that, 
it makes it easy to sum up the bubble expansion of the 
p-h residual interaction (for a review see Ref. 0 )- An ad¬ 
ditional advantage is that one is incorporating the short- 
range p-p correlations in the p-h effective interaction. In 
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terms of the Landau parameters one gets 




F 


G 


l + A{q)F l + A{q)G 

G' 


[cT ■ cr') 


1 + A{q)G^ 


-{a ■ -t’), (4) 


where F = N{0)f, G = N{0)g, F' = N{0)f', G' = 
N{0)g' and N{0) is the level density on the Fermi sur¬ 
face. A(g) is the Lindhard function resulting from the 
integration of the polarization insertion [iSllsee Appen¬ 
dices). This interaction will be adopted as vertex inser¬ 
tion (wiggly line) in the two-bubble diagram as shown 
in Fig. So it represents the missing RPA correction 
in the one-bubble approximation for the screening, al¬ 
ready calculated with the Gogny force 0]. Of course, 
this choice entails that the Landau parameters must be 
calculated using the same interaction, i.e. the Gogny 
force in our calculations, adopted for the two external 
interaction vertices. The Landau parameters with the 
Gogny force D1 |^, after RPA bubble summation at a 
given density, depend on energy and momentum transfer 
via the Lindhard functions and their asymptotic values 
coincide with the previous values at the same density. 

At this point let us again discuss our choice of the 
Gogny force in the vertices (dashed lines) of Fig. Q] In 
principle a good approximation to the dashed lines in 
Fig. lb and Ic would be a microscopic G-matrix. This is 
well known from many body theory. We here replace the 
G-matrix by the phenomenological Gogny force which 
has precisely been adjusted to a G-matrix calculation (see 
Ref. ). More questionable is our use of the Gogny force 
in the Born term of Fig. la. In principle in Fig. la the 
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bare force should be taken and the use of the Gogny force 
is for pure convenience here because a bare force scatters 
to very high energies. However, for better quantitative 
comparison of the Born terms, Fig. la, and the screening 
terms. Fig. lb and Ic, rather the use of an effective low 
energy force like uiow-k 0, which is phase shift equivalent 
with the bare force at low energies, would be appropriate. 
Since for us such a renormalised uiow-kfs not available in 
this work and because we know from Q that uiow-k and 
Gogny forces are quite similar, we feel entitled to use the 
Gogny force for the Born term of Fig. la. In view of the 
fact that, as will be shown below, the contributions of 
the screening terms are not at all small compared with 
the Born term, we feel that a slight inaccuracy in the 
evaluation of the Born term will not at all have strong 
consequences for the conclusions which will be drawn in 
this paper. In this respect we also should mention that 
in the channels, the density dependent part of the 
Gogny force drops out in the Born term. However, this 
is not the case in the ^Si channel. 

Different screening mechanisms come into play accord¬ 
ing to whether the medium is nuclear matter or neutron 
matter, and also whether one has pairing between like 
or unlike particles. Denoting by S and T total spin and 
total isospin in the p-p channel, the diagrams of Fig. ^ 
are written in the spin and isospin representation 

^ST = VU + Vf^ + Vi% (5) 


where 


(TiTi 


= X! (-l)'^^’'''^^C'(CT-CTi;cr'CT2|5'i)G(T-Ti;T'T2|Ti) X G(tT2a-;-0-i0-'|*S'2)G(T2T;-TiT'|r2) 

(JiTiSiTiki 

X {kki\Vsj^Ti\k'k2) {k2k\Vs^T2\kik') A{ki,k2)GsT, 

= X! - cri;cr'cr2|5'i)G(r - ri;r'T 2 |ri) X G(cr4d;-cr3d'|5'3)G(T4f;-r3f'|r3) 

(TiTiSiTiki 

X {kki\Vs,T,\k'k2) Vifrt {k4MVs,n\kJ^') A(fci, fe)A(fc3, A:4 )Gst, 


in which 

Gst = C{aa-a'a'\S)G{Tf-T'f'\T) 

and 

G(tTiCr2;cr3fT4|5') = S Ui + (J 2 ) 

X ^2^32^4] S (T3 -t- (J4) (jcTi-t-cr2,cr3-t-cr4 ■ 

It is worth noticing that the insertion is depending 

only on the total spin and isospin in the p-h channel, 


whereas the two external interactions are expressed in 
the spin and isospin coupling mixed representation of p- 
p and p-h channels. 

Nuclear matter within the RPA treatment of p-h resid¬ 
ual interaction suffers from a mechanical instability sig¬ 
naled by the well known singularity of the compression 
modulus at Fq = — 1. It occurs in a density domain 
where pairing gap is large, and its influence on the screen¬ 
ing compromises any quantitative prediction. In neutron 
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FIG. 2: Splitting of the ^Sq pairing screening potential Vkpk'i^) in nuclear matter. Left (right) panel is for the neutron (proton) 
p-h one-bubble insertion. The Fermi momentum is fixed to 1.0 fm“^. 


matter the instability would not appear, but for several 
interactions including Gogny Dl, Fq approaches danger¬ 
ously — 1 so making doubtful any estimate. This draw¬ 
back has been cured by the Babu-Brown theory 
where the class of bubble diagrams is renormalized using 
the concept of induced interaction. The latter is implic¬ 
itly dehned as follows 

hp/i — krfir T ^ind — ^dir “t” '^KPAiy^ph'): (6) 

where the first term Vdir is the direct residual interaction 
(G- matrix or Gogny force in the simplest approximation 
[l9|) and the second one is the RPA bubble series where, 
the vertex insertions are given by Vph itself instead of 
Vdir- The construction of the induced interaction ap¬ 
proximation (IIA ) is, in general, a very complex problem 
[ll Il9ll^ 1^ 1^ . but it is easily realized in terms of Lan¬ 
dau parameters, as shown in detail in the Appendix E. 
In principle, the propagators and the RPA calculations 
should be calculated in the superfluid matter. However 
we may guess that this consideration induces only a sec¬ 
ond order effect and we study the pairing screening only 
in the normal system in the present paper. 

III. RESULTS 

Based on the approximations discussed in the preced¬ 
ing section, the screening interaction has been calculated 
for the nn (or pp) isospin-triplet channel and np 
isospin-singlet channel. In the former case the screen¬ 
ing effects due to both pure neutron matter and sym¬ 
metric nuclear matter were also estimated for the sake of 
comparison. In the past the study of isospin-triplet pair¬ 
ing in nuclear matter has usually not been considered, 
since one mainly had in mind the problem of superflu¬ 
idity in neutron stars, but for applications to nuclei the 
investigation of pairing in nuclear matter is also relevant. 

As in a previous study Q , we used for the calculations 
the Gogny D1 force (see Ref. and Appendix A for 


details) for the external vertices, as well as for the Landau 
parameters describing the p-h residual interaction in the 
internal vertices. 


A. Screening potential in the isospin triplet 
channel 

Let us first consider the screening of symmetric nuclear 
matter and neutron matter on pairing between two nucle¬ 
ons in the channel. In the one-bubble limit the neu¬ 
tron matter polarization gives rise to a repulsive screen¬ 
ing of the nn channel, as shown in Fig. 3 of Ref.|^. 
This is a well known result since long , and it 

has been interpreted as due to the dominance of the spin 
density fluctuations over the density fluctuations . It 
results in a sizeable quenching of t he p airing gap, rein¬ 
forced by the self-energy effects [3,123, lla- In nuclear 
matter, the screening of nn pairing in the channel 
can be split into two parts: neutron bubble insertion and 
proton bubble insertion, as shown in Fig. |21 The neu¬ 
tron bubble produces the same effect as for nn pairing in 
neutron matter, but here the additional proton polariza¬ 
tion gives a very large attractive contribution, as shown 
in the right side of Fig. Inland in Fig. |21 As a result the 
full medium polarization is enhancing the isospin triplet 
pairing (anti-screening), and produces, at least in the 
one-bubble approximation, a huge value of the pairing 
gap [13, which is cancelled only partially by the disper¬ 
sive effects of the self-energy treated on the same footing. 

In Fig.|21we show the full w dependence of the screen¬ 
ing part of the NN interaction for the intermediate Fermi 
momentum kp = 1.0 fm“^ and in Fig.|3we show the di¬ 
agonal part of the sum of the neutron and proton bubble 
insertions for k = k' = kp with kp ranging from 0.6 to 
1.4 fm“^. The total (i.e. Born plus screening on shell nn 
pairing interaction on the level of one bubble exchange 
is shown in Fig. for the channels in nuclear and 
neutron matter as a function oi kp (dashed lines). In 
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FIG. 3: Diagonal part of the sum of the neutron and proton 
bubble insertions in nuclear matter with kp ranging from 0.6 
to 1.4 fm“^ 


spite of the slight decrease of the intensity of the induced 
interaction at very low nuclear matter densities {kp < 
0.5 fm“^), it still seems to considerably enhance the at¬ 
traction of the Born term (solid line) in the low density 
limit. This qualitatively seems to be in line with the find¬ 
ing of Heiselberg et al. where it was found that for a 
four component Fermi system (nuclear matter) the gap 
should be enhanced in the zero density limit, whereas, 
on the contrary, in a two component Fermi system (neu¬ 
tron matter) the gap should be reduced. This is true if 
the scattering lengths in all channels are approximately 
equal. This is not the case in nuclear matter where the 
T = 0 scattering length is a factor three to four smaller 
that the T = \ one and therefore a more accurate inves¬ 
tigation of the range 0 < kp < 0.5 fm“^ is still in order 
but numerically very delicate. 



hioiMeV) 


FIG. 4: Nuclear matter screening potential Vkpy{^) in the 
channel in the one bubble approximation. The Fermi 
momentum kp is hxed to 1.0 fm“^ 


B. Screening potential in the isospin singlet 
channel 

In nuclear matter one also has to consider the isospin- 
singlet pairing, namely the pairing between unlike nu¬ 
cleons. The dominant coupling interaction is due to the 
^SDi component of the force (in the following we ne¬ 
glect the D component, which makes the calculations 
very complicate giving however a small effect). As men¬ 
tioned above, the magnitude of the gap in this case is, 
using the bare force, too large with respect to the values 
observed or predicted in finite nuclei [3, 111) ■ Then 

one may expect that the screening is reducing the gap 
to a more physical value. And in fact it turns out to 
be repulsive, at least in the one-bubble approximation, 
as shown in Fig. 0 Actually there are energy domains 
where the screening potential is attractive, but it is re¬ 
pulsive around tu = 0, which is most relevant for pair 
formation. This repulsive effect is also confirmed for the 
full on shell pairing force in Fig. 0 especially towards 
lower densities (dashed line on left panel). 

C. Screening with induced interaction 

Regardless of the physical implications of the screen¬ 
ing effects on the pairing, what we learn from the above 
results is that the order of magnitude of the screening 
or anti-screening is not small at all. This indicates that 
the full bubble series (RPA) should be summed up. On 
the other hand, we know that a pure RPA approxima¬ 
tion, namely the bubble series with the Gogny force as 
residual p-h interaction, is singular in the instability re¬ 
gion. The only way to remove this drawback is to intro¬ 
duce the induced interaction. So far the full RPA with 
induced interaction has been applied to neutron matter 
using a G-matrix instead of the Gogny force P| . For the 
present calculations we adopted the same approximation 
except that, as discussed in the preceding section, the 
p-h residual interaction is expressed in terms of Landau 
parameters. This makes the resummation of the bubble 
series much easier to perform. 

1. Landau parameters 

According to the motivations of the preceding section, 
the screening interaction due to the p-h excitations is cal¬ 
culated approximating the residual interaction with the 
I = 0 Landau parameters (see Ref. and Appendix B). 
The results are plotted in Fig. 0 both for nuclear mat¬ 
ter and neutron matter. As the Brueckner G-matrix, the 
Gogny force, where the short-range correlations are also 
incorporated, does not prevent the mechanical instability 
{Fq < — 1) to appear in the low density domain of nuclear 
matter close up to the saturation point. Otherwise the 
values of the Landau parameters with the Gogny force 
reproduce satisfactorily the empirical values. In neutron 
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FIG. 5: Landau parameters from the Gogny force Dl. Left 
panel is for nuclear matter and the right one for neutron mat¬ 
ter. 






FIG. 7: Energy(ftoi) and momentum-transfer(g) dependence 
of the induced Landau parameters (see Appendix E: E2). The 
upper two panels are for neutron matter and the lower four 
for nuclear matter. The Fermi momentum kp is fixed to 1.0 
fm-^ 


imation, whereas asymptotically \tujj\ 3> 0 they tend to 
the zero order limit, shown in Fig. El The dependence 
on the momentum transfer is not significant in the range 
shown in the figure, but in a wider range it could. 


FIG. 6: The same as in Fig. El but in the induced interac¬ 
tion approximation. Only the values at the Fermi surface are 
depicted. 


matter Fq slowly decreases, since the Gogny Dl force 
misses the repulsive rearrangement term in the neutron 
channel. The singularity due to the mechanical instabil- 
ity is removed by the Babu-Brown induced interaction 
as shown by the enhancement of Fq in Fig. El Since 
the induced interaction entails a strong coupling among 
the components of the p-h interaction in nuclear matter 
(see Appendix D), the rapid rise of Fq makes the other 
Landau parameters to rapidly bend up also. Despite the 
simplicity of our approximation the induced interaction 
still keeps dynamical effects because of its dependence 
on the p-h propagator. In Fig. [Tithe Landau parameters 
versus energy and momentum transfer in the induced in¬ 
teraction approximation, as shown in Appendix E, are 
plotted. The Huj = 0 and q = 0 values coincide with the 
static values of Fig. El in the induced interaction approx- 


2. Screening potential from the full RPA residual 
interaction 

Using the Landau parameters as vertex insertions in 
the p-h bubble diagrams the screening potential has been 
calculated in the full RPA limit with induced interaction. 
The pure RPA, i.e. without induced interaction, has not 
been considered, since it displays a singular behavior in 
nuclear matter. Fig. El summarizes the main properties of 
the full pairing interaction in comparison with the Gogny 
force. In the isospin triplet channel (^^o) in neutron 
matter the net screening potential is repulsive because of 
the dominance of the repulsive spin fluctuations over the 
attractive density ones. A suppression is in accordance 
with practically all existing calculations, see e.g. |^. 

For instance we also see from Fig. El that a substan¬ 
tial suppression survives at low densities. This is in 
agreement with the recent studies of Gori et al.ls^ and 
Schwenk et al. M where also an important reduction of 
the gap in low density neutron matter was found. Gom- 
paring with the low density behaviour of the one-bubble 
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FIG. 8: Diagonal matrix elements {k — k' = kp) of the interaction at hu) = 0: Born term with Gogny force (solid lines), 
Born term with one bubble insertion (dashed lines) and fnll pairing potential including one and higher bnbble insertion (IIA) 
(dash-dotted line). The baryon density corresponds to kp = l.Q fm“^. The Born term with the Gogny force (solid line) in the 
^So channel is density independent, whereas it is density-dependent in the ^Si channel. 


induced interaction (see Fig.4 of Ref. [3]), we learn now 
that the one-bubble approximation is insufficient. At 
kp = 0.3 fm“^ its contribution to the pairing potential is 
7.93 MeV, while the bare interaction is still -36.7 MeV. A 
sizeable effect is expected in the limit N{0)V <C 1, which 
for the nuclear force in the T = 1 channel (scattering 
length -18.8 fm) is fulfilled only at densities much less 
than 10“^ fm“^, i.e., much beyond the range of interest 
of the present work. 

In nuclear matter the dominance of the attractive force 
exerted by the proton environment on nn pairs (or neu¬ 
tron on pp pairs) which we already found at the one bub¬ 
ble level, becomes more pronounced in the full screening 
potential. This dominance becomes dramatic in the case 
of the ^5'o channel. The potential plays a role of strong 
anti-screening indeed. But its effect on the pairing corre¬ 
lations should be counterbalanced by the self-energy ef¬ 
fect which in nuclear matter are much more pronounced 
than in neutron matter. The amount of compensation 
can be seen in Table 1. In the case of spin triplet pairing 
there is a sizeable compensation between self-energy and 
vertex corrections at any density. For the spin singlet 
pairing in neutron matter, as expected, both quench pair 
correlations, except around kp = 1.0 fm“^. The self¬ 
energy is also reducing pairing in nuclear matter, 
however, the T = 0 part of the induced force is enhanc¬ 
ing the pairing interaction very much, so that a strong 
(unphysical?) increase of the gap can be expected. 

IV. DISCUSSION AND CONCLUSION 

In this study we included medium polarization effects 
in addition to the Born term (Fig. la) in infinite matter 
in three different channels: i) pairing in pure neutron 
matter, ii) pairing in symmetric nuclear matter, iii) 
^SDi pairing in symmetric nuclear matter. As already 


TABLE I: Comparison of dimensionless pairing interactions 
in three approximations: bare interaction (first column), self- 
energy effects (second column), self-energy plus screening ef¬ 
fects (third column), “np” in the second row and “nn” in 
the third row stand for nuclear matter and neutron matter, 
respectively. The factor Z is dehned in Eq.(9) of Ref. @1. 


channel 

fcF(fm 


Ao‘ZVbare 

+ Ucr.) 


0.7 

-1.21 

-1.04 

-1.47 


1.0 

-0.95 

-0.54 

-0.70 


1.3 

-0.21 

-0.15 

-0.23 


0.7 

-0.43 

-0.38 

-2.46 

^S'o(np) 

1.0 

-0.41 

-0.23 

-1.05 


1.3 

-0.31 

-0.22 

-0.33 


0.7 

-0.44 

-0.33 

-0.26 

^S'o(nn) 

1.0 

-0.41 

-0.33 

-0.26 


1.3 

-0.31 

-0.26 

-0.21 


outlined in the introduction, one expects from the polar¬ 
ization contributions (self-energy -|- vertex) quenching of 
pairing in i), slight additional attraction in case ii), and 
quite strong quenching in iii). From Fig. |Hlwe can see 
in how much these expectation are realized within our 
calculations. We investigated two approximations. The 
first just consists in considering the one bubble exchange 
between two nucleons (Fig. lb). The second also takes 
into account interactions between the particle and the 
hole (RPA). These latter interactions, are approximated 
by Landau parameters, renormalised by the Babu-Brown 
procedure (Fig. Ic). 

The results in neutron matter for ^S'o pairing confirm 
the ones of other authors, i.e., the bare interaction is 
screened by 20 to 30%. This screening is stronger in the 
one bubble approximation than in RPA. As expected, 
medium polarization gives extra attraction in the ^S'o 













channel in symmetric nuclear matter (case ii)). However 
the extra attraction turns out to be enormous. With re¬ 
spect to the one bubble approximation the resummation 
of the bubbles in RPA, still enhances the effect. With 
such a strong total pairing force the gaps in finite nu¬ 
clei would clearly be out of scale, when applying an LDA 
scenario. One is wondering what is happening. Maybe 
LDA is quantitatively completely wrong when going from 
nuclear matter to finite nuclei. In the past we have 
shown 1^ that for RPA, LDA is quite good for fi¬ 
nite momentum transfers. Here the momentum transfers 
are concentrated around <7 = 0 where LDA is not quite 
valid. This may be the reason. A second uncertainty 
comes from the theoretical approach. We see from Fig.|HI 
that going from one bubble to RPA, the effect is strong 
and therefore the many body approach may not be con¬ 
verged. Though we very carefully checked our formulas 
and the numerical procedure, an independent confirma¬ 
tion of this surprisingly large effect would be indicated. 
In any case it hints to the fact that also in finite nu¬ 
clei anti-screening should be taken seriously as this has 
indeed been the case recently 

In the ^Si channel of nuclear matter we expect quench¬ 
ing of the pairing force due to medium effects. This is 
indeed what is happening in the one bubble approxima¬ 
tion. However, as in all other channels, additional bubble 
summation yields extra attraction, so that with RPA also 
in the channel one obtains slight extra attraction. 

So let us summarize the situation. At least in the one 
bubble approximation the trends in all three cases are 
as expected: repulsion for cases i) and iii), attraction for 
case ii). In all these cases adding RPA bubble resum¬ 
mation with respect to the one bubble approximation 
leads to considerably more attraction. In case iii) this 
reverses the expected trend from additional repulsion to 
additional (slight) attraction (with respect to the Born 
term). Whether these additional RPA correlations go 
into the right direction is not clear. Their strong effect 
means that the many body scheme is not converged. An 
approach based on a more controlled scheme like a vari¬ 
ational theory would be in order. Also the question in 
how much the present study can be linked to finite nu¬ 
clei needs further studies. We, however, believe that at 
least the trends should be the same in nuclear matter and 
finite nuclei. 

In the end let us again comment on our use of the 
Gogny force in the Born term. As already mentioned 
this is not completely consistent since for densities larger 
than « po /4 the Gogny force yields a considerably larger 
gap than the bare NN force (which should be used in full 
rigor). However, in this work we are only interested in 
relative trends of the Born term versus the induced force. 
Since the latter turns out to be a sizeable fraction of the 
Born term, a small change of the lowest order term will 
not invalidate our conclusions of this work. On the other 
hand we feel that the situation is not sufficiently under 
control to warrant the huge amount of numerical work 
to calculate the gap in the different channels. While in 


neutron matter this may yield reasonable results which 
we will give in a future publication, we can already say 
that in the other channels and for instance for T = 1 
in nuclear matter the gap values will be extremely large 
invalidating the whole weak coupling approach of BCS. 
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APPENDIX A: GOGNY FORCE 

The appendix is devoted to develop the formalism de¬ 
scribing the interaction and the approximations adopted 
in the calculations. The Gogw force, neglecting the spin- 
orbit term, is given by Ref. [g 

{hk2\V\k3k^) = ^ 

i^l,2 

-H,P^ - M,P^P^) 

+7(1+ (Al) 

where k = (fc, a, r) denotes the single-particle state and 
/i, 7 is defined as 

7 = (A2) 

which will be used in Appendix B, C, D. The parameters 
of Gogny D1 force are reported in Table 1. 


TABLE II: The parameters D1 for Gogny force, 
(ts = 1350 MeV fm^) 



ri (fm) 

m 

Bi 

Hi 

Mi (MeV) 

1 

0.7 

-402.4 

-100 

-496.2 

-23.56 

2 

1.2 

-21.3 

-11.77 

37.27 

-68.81 


APPENDIX B: LANDAU PARAMETERS FROM 
GOGNY FORCE 

The p-h residual interaction can be described in terms 
of the Landau parameters to be extracted from the pre¬ 
ceding Gogny force. In order to do that we introduce the 
parameters 

nka,T (Bl) 
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where 

{ 1 baryon density (5=0,r=0), 

(—spin density (5=1,T=0), 

(—isospin density (5=0,T=1), 
(—1)'^+'^“^ spin-isospin density (5=1,T=l). 

Spin and isospin asymmetric nuclear matter can be is 
expressed in terms of the previous parameters as follows 

n=lC + ]-^{kk'\V\kk')n',,nk', (B2) 

^ kk' 

where K, is the kinetic part and V denotes anti¬ 
symmetrized interaction. Then the Landau parameters 
are obtained as functional derivatives 

Fi^' = = iV(0) ■ (B3) 

The normalization factor N(0) is the level density and 
it is introduced to make the Landau parameters to be 
dimensionless. After expanding into Legendre polynomi¬ 
als, we select for our purpose only the zero-order Landau 
parameters. For nuclear matter we get 

i 

-2H, - 4M,)u,] + |7, 

i 

Go = ^[(2B,-M,)«.-(VF,-2iJ,Vz] + H 

i 

G'o = ^[-M,u,-WiU,]-i7, (B4) 

i 

where Ui = ^iN(0)/A and vt = Uie~^'(sinh. Zi)/Zi and 
Zi = {rikpY/2 and |fc| = |fc'| = kp, kp being the Fermi 
momentum. For nuclear matter iV(0) = 2m*kp/{pfif'. 
7 and /i^ are defined in Eq. (A2). 

In the case of pure neutron matter the zero range term 
of the Gogny force does not contribute. Only two (zero 
order) Landau parameters survive, i.e. 

Fo = ^[(2W', + F,-2iJ,-M,)uz 

i 

-{W^ + 2B,- H,-2M,)v,], 

Go = J2[{B, - M,)u, - (W, - H,)vi], (B5) 

i 

where now Ui=tJLiN{0)/2 and N{0)=m*kp/{phY is the 
level density of neutron matter, kp being the neutron- 
matter Fermi momentum. 


APPENDIX C: ONE-BUBBLE SCREENING 
INTERACTION 

The medium polarization effects on the NN interaction 
in the RPA limit are given by the bubble series. The one- 
bubble term (diagram) has the following expression, 

= -Y, {kki\V\k'k2) {k2k\V\kik') 

kik2 

(Cl) 

which is suitable for the pairing interaction where one 
particle in the state k = {k,a,T) couples to one parti¬ 
cle in the state k = {—k, a', t'). The function L is the 
polarization part M 

Lkik2 {^) = J {^')Gk2 (C 2 ) 

After w-integration within the single-pole approximation 
for the Green’s function, one obtains 

j / \ 7rfc,^(l nk2)Zki 

Fkik2{^) = -;-;-- H-;- 

t^/ci + u;k2 + ^ - ^ki + - w - ZT? 

where the Z-factors embody the self-energy effects. In 
the coupled case, Eq. m becomes, 

V^l\k,k',uj) = -^ JY (ST)Lfe,fc,(a;)dL!fcdIIfcq 

kik2,ij 

(C3) 

where ST stands for ^5o or ^5i channel, and z, j takes 1 
and 2. In the neutron matter medium, 

v}pCSo) 

= 2(a, - l3i){B, -H,-M, + VF,)[/3,(F, - IF,) 
+7,(5, - M,)] + [a,{H, - IF,) + /3,(5, - M,)] 

X [/?, (M, - 5,) + 11, (IF, - 5,)], (C4) 

where Q!,=^iexp[-irf (fc - fc 2 )^], /3i=/Ziexp[-irf (fc - 
^0^]; ^i=+iexp[—irf (A; — fci)^], and in nuclear matter 
medium, 

C,fC5o) 

= —{ciiHi + fiiBi + ^^)(PjBj + rjjHj + ^ 7 ) 

+2(a, - A)(5, + IF,)[/3,(5, - IF,) 

+7,(5, - M,)] + 2[a,(5, + M,) + /3,(5, + IF,) 
+'y][PjWj + iijMj + ^ 7 ] - [a,(IF, - Hi) 

+/3,(M, - 5,)][/3,(M, - 5,) + VjiW, - 5,)]. (C5) 

In ^5i channel, 

v^,esi) 

= -{aiB^ + /3,5i + ^j){/3jH, + 7 , 5 ,- + ^ 7 ) 

+2[a,(M, - 5,) + /3,(IF, - 5,)](/3,IF, + 7 ,M, 

+ ^ 7 ) + 2(a, - /3,)(M, - IF, + 5, - 5,)[/3,(5, 

+IF,) + r]j(Hj + Mj) + 7 ] — [ai{Bi + IF^) + Pi{Hi 
+Mi) + 7 ] [j3j {Hj + M,) + 7 , (5, + IF,-) + 7 ], (C 6 ) 

where 7 is density dependent, as defined in Eq. (EH). 
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APPENDIX D: RPA SCREENING 
INTERACTION 

The summation of the full RPA bubble series can be 
carried out only in few cases as explained in Sec. II. 
A reasonable approximation which preserves the correct 
evaluation of vertex insertions between p-p external lines 
and the p-h internal lines is depicted in Fig. 1(c). It is 
based on a two-bubble approximation where the internal 
vertex is replaced by the full RPA p-h interaction 

p(3) = ^ {kki\V\k'k2)Lk,kA^) 

kik2ksk4 

x{k2k3\VB.pA\kik4:)Lk3ki{^) {ki^V\k3k') . 

(Dl) 

The momentum transfer in the pairing interaction is 
q = k — k' = k 2 — ki = k^ — k 3 . Using the Landau 
parameters, Urpa is only depending on q. The latter 
is calculated using the Landau parameters for which the 
RPA summation can be easily performed. In turn, this is 
a reasonable approximation, since the relevant p-h exci¬ 
tations are those nearby the Fermi surface where the Lan¬ 
dau parameters are calculated. Therefore, in the Landau 
approximation, 

^RPA = f + gO-l-(T2 + f'Tl-T2+g'{cTl-(T2){Tl-T2), (D2) 

where /, g, /', g', including so called induced interaction, 
will be given in Appendix E. In the coupled case, the 
interaction with full RPA (the bubble series summation 
with two and higher bubbles) is 

^ j dVLkdVLk- Yl (D3) 

k2k4,ij 

where ST could be ^Sq or 35*;^. For channel, Uy takes 
the form 

+ Q^|di){P]aj + Qjqj)f 
-3(X,a, + + Yjf]j)g (D4) 

in neutron matter, where at = /iiexp[—r^(A; — A;')^/4], 
Pi = /riexp[-rf(fc-^ 2 )^/ 4 ], qt = fiiexp[-rf{k + k4)'^/4], 
P^ = B,- 2H, -M, + 2IU„ Q, = -2B, + iL, + 2M, - IU„ 
Xi = B^- M^, Y^ = H,- Wi, and 

V^P = {Aiar - EiPi + l'^){Ajaj - Ej^j + ^j){f - f) 
-3(Cjaj -H A/3 z + - g') 

(D5) 


in nuclear matter where Ai = 2Bi — 2Hi — Mi+AWi, Ei = 
2B, - 2H, - 4M, + W,, Q = 2B, - A = 2A - W,- 
For 3 A channel, we have 

Vlp = (Aiai - EiPi Y ^'Y){Ajaj - E^q^ + |7)(/- /') 

+ (Ciai + DiPi + + Djqj + ^^){g — g')- 

(D6) 

APPENDIX E: INDUCED INTERACTION 


The induced interaction is given by Eq. ®, i.e. 

^^RPA = Vdir + Vind- (El) 


where the induced interaction Vind is the bubble series 
with the full interaction itself inserted in the interaction 
vertices. Its explicit expression, in nuclear matter case, 
is. 


f — fd + find, f — fd + //nd> 

g = gd + 5ind, g' = g'dY 1 ?^: (E2) 

where 


4/ind — Af + 3Ag + 3Afi + 9Ag', 

4<7ind — Aj Ag Y 3Ay/ 3A^/, 

4/ind = ^/ + 3Ag — A// — 3Ag/, 

45i'nd = Af-Ag-Af,YAg,, (E3) 


and 


A„ 


A 


1 + A(0)/A 

Hq) ~2 


/^ ^/' = 


A 


1 -t iV(0)/'A 




1 -h N{0)gA 


g I “ 


A(<Z) 


1 -f N{0)g'A 


jrr- (E4) 


While in neutron matter case, only / and g is remained 
in Eq. and the self-contained equation is 


2/ind — Af Y 3Ag, 

25ind = Af - Ag, (E5) 

where A / an d Ag has the same form as in Eg. (IE3II . Usu¬ 
ally En. dEAll and IIE5II is solved iteratively, in which the 
first guess is the direct term fd- The simplest approxi¬ 
mation to calculate the induced interaction is given ex¬ 
pressing / in terms the L = 0 Landau parameters where 
A = A{q, 0) is the Lindhard function in the static limit 
a; = 0. 
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